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We study a family of closed quantum graphs described by one singular vertex of order n = 4. By 
suitable choice of the parameters specifying the singular vertex, we can construct a closed sequence 
of paths in the parameter space that physically corresponds to the smooth interpolation of different 
topologies - a ring, separate two lines, separate two rings, two rings with a contact point. We find 
that the spectrum of a quantum particle on this family of graphs shows the quantum anholonomy. 

PACS numbers: 03.65.-w: quantum grapli, boundary condition, eigenvalue anholonomy 



INTRODUCTION 

The idea of quantum mechanics on graphs has been in- 
troduced in the last century. During the last 30 years, it 
has been developed into a powerful tool for the study of 
particles on tiny graph-like structures [5]. Since a graph 
is a one-dimensional variety, the analysis of a particle on 
a graph consists in solving a system of ordinary differ- 
ential equations, which makes the graph models simple 
from the mathematical point of view. Therefore, quan- 
tum mechanics on graphs can serve also as a useful lab- 
oratory for the study of various complex quantum phe- 
nomena in a simple setting. 

Here we examine the influence of the topology of quan- 
tum graphs on spectral properties, taking inspiration 
from the works [T] and [7] . We may regard that the topol- 
ogy change mimics a macroscopic or violent variation of 
systems, e.g., the rupture and connection of quantum 
wires. We will show that a parametric evolution of an 
eigenenergy may form an open trajectory along a cycle 
involving the topology changes of graphs. If this is the 
case, the trajectory of the corresponding eigenspace also 
forms an open cycle. Such a discrepancy of eigenobjects 
is called quantum anholonomy [H IH [H] • 



THE MODEL 



sticked together at their outer endpoints to form a ring 
of length L2, giving the whole object the shape of the 
character "00" (Fig. [l]). We assume there is no poten- 
tial on the lines outside the central vertex. Therefore, 
a quantum particle confined to this object moves freely 
on the lines, and the only non-trivial part of the system 
is the connection condition of the wave function at the 
node. 

We assign the coordinates xi e [0, Li] and X2 € [0, L2] 
to the left and the right ring. The values xi = and 
X2 = correspond to the inner endpoints 1 and 3, re- 
spectively, whereas the values Li and L2 correspond to 
the inner endpoints 2 and 4, see Fig. [T] The behavior of 
a particle in the central vertex of degree 4 is determined 
by a boundary condition in the vertex. The general form 
of a boundary condition in a vertex of degree n follows 
from the study of self-adjoint extensions of the Laplacian 
operator on a graph. The boundary condition consists of 
n equations connecting the boundary values and deriva- 
tives that can be written in the form [6] 



+ B*' = 0, 



(1) 



where ^P, ^' are vectors formed from the values of the 
wave function in the vertex and from its derivatives 
(taken in the outgoing sense), respectively, and A and 
B are n x n complex matrices with the property 




FIG. 1. A closed quantum graph with an n = 4 singular 
vertex. 

Let us consider a singular node with four outgoing one- 
dimensional lines of finite lengths. Let the first two lines 
be sticked together at their outer endpoints to form a 
ring of length Li, and likewise, the second two lines be 



AB'' = BA"^ and T:ank{A\B) = 



(2) 



with {A\B) denoting the n x 2n matrix composed from 
the columns of the matrices A and B. For the graph 
depicted in Figure [l] we have 



/ V'i(O) \ 

MO) 



^'i(O) 
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\-ML2)J 



(3) 



In our model, we focus on a specific subset of all possible 
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connection conditions ([T]), namely 



/I tb{e) 

l-c{e) 

\0 



ta{e) sd{0)\ 

ta{e) 

c{9) tb'(0) 

/ 

/ 0^ 

-tb{e) c{0) 

-ia(6') l-c{0) 



* (4) 



where t and s are positive constants, and a{9), b{9), b'{d), 
c{6) and d{9) are real periodic functions of 9 with period 
27r that vary within the range between and 1. In ad- 
dition, it can be easily shown that the functions must 
satisfy 



a{9) 7^ 
b{9),b'{9)^Q 



c{9) = 0, 
cs{9) = 1 



(5) 
(6) 



for all 9 e [0, 27r] in order that the matrices A, B involved 
in the connection condition Q obey 



Long Cycle 

At first we consider the boundary conditions Q with 
fmictions a{9), b{0), b'{9), c{0) and d{9) chosen as 



a{9) — cos ( 

m 

b'{9) = (sine - V3COS0 - V3 



cos 

2 



cos 6* + - 



c{9) 




sin 9 — \/3 cos 9 — \/3 | , 

• ^1 1 

sm 

2' 2 



. V3 

sm — 

2' 2 



(7) 



(8) 



d{9) = ^ ||sin6l| -sin6' + \/3 



sin 9 - 



V3|} 



(9) 
(10) 



The graphs of these functions for € [0, 27r] are plotted in 
Fig. [2] The interval [0, 27r] can be divided into six regions, 
in each of them the vector {a{9), b{9), b'{9), c{9),d{9)) has 
a different structure of zero and nonzero entries. As we 
show below, these regions correspond to mutually differ- 
ent topologies of the system. It extends the basic idea of 
parametric systems with a varying topology introduced 
in works of Balachandran et. al. [T] and also in [7] . Our 
use of a single n = A singular vertex makes the treatment 
more systematic and unified. 



VI V VI 
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FIG. 2. Graphs of functions defined by Q for £» £ [0, 2-k] 



Region I For 9 £ [0, we have 



/I t-a{9) \ 
1 t ■ a{9) 



't-a{9) 1 

-t-ai9) Ij 



vf, (11) 



where ta{9) takes the value i at 6* = and it continuously 
decreases down to at = ^ . With regard to the struc- 



ture of (11), the endpoint 1 is joined to the endpoint 2 



and the endpoint 3 is joined to the edpoint 4. On the 
other hand, the endpoints 1 and 2 are disjoint, so are 3 
and 4. Therefore, the system is topologically equivalent 
to a single ring, see Fig. |3] 

Region II For 9 e [|, ^], we have 



/I 



1 - c{9) 



c{9) 



c{9) 



1 - c{9) 



\ 



(12) 



V 



where we start with the value c(^) = and continuously 
increase it up to c(^) — 1. Since all four inner endpoints 
are disjoint, the system is topologically equivalent to two 
separate lines. 
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Region III For 6 e [^^^jTt], we have 



Region VI For 6* e 



we have 



/I t-b{e) 



\ 



1 t-b'{e) 
( 



-t-b{e) 1 



V 



(13) 



-t-b'{e) ij 



to tb{-K) = t, while i5'(6') takes the value tb'{0) = at 
6 — continuously increases up to i at 6* = ^ and 
conies down back to at = tt. The connection condi- 



tion (13) implies disjoints endpoints 1 and 3, as well as 
2 and 4, making the system topologically equivalent to 
two separate rings. 

Region IV For 6 € [n, 4^], we have 



/l t-b{e) s-d{9)\ 



-t-b{e) 1 

V-s • d{9) l) 



vf, (14) 



where tb{6) continuously decreases from tb{^) = t to 
tb{^) — 0, while sd{0) continuously increases from 
sd{^) = up to sd{^) — s. The system is topologi- 
cally equivalent to a ring and a line connected at a single 
point since three of four endpoints, namely 1, 2 and 4, are 
connected in the node. Note that at the juncture of re- 
gions /// and IV, the system is topologically equivalent 
to a separated line and a ring. 

Region V For G we have 



/I t-a{0) s-d{d)\ 
1 t ■ a{0) 



V 



-t ■ a{9) 1 
\-s-d{9) ~t-a{0) Ij 



(16) 



where tb{0) continuously increases from tb{^) = up ta{2n) 



where ta{0) continuously increases from ta(^) up to 



sdif) 



t, while sd{0) continuously decreases from 
s down to sd{2TT) — 0. In this region, the 
system is topologically equivalent to two rings connected 
at a single point since all four endpoints are connected 
together at the central node. 



V 




FIG. 3. Cycle of topological changes corresponding to func- 
tions Q. 

Therefore, with the parameter change = 27r, 
seven distinct topologies are traversed in a continuous 
manner. The situation is illustrated in Fig. [3] The sys- 
tem is never identical for any two different values of 
in [0, 27r), because the the functions a{0) to d{0) form 
district combinations of patterns in all six regions. 



/I 



1 ' c{0) 



c{0) 



( \ 

c{0) 

1 - c{0) 

Ij 



(15) 



where c continuously decreases from c(- 



1 down to 



3 ' 

c(f^) = 0. The syste m is topologically equivalent to a 
single line, since only endpoints 1 and 4 are connected, 
whereas 2 and 3 are disjoint. 



Short Cycle 

As an additional model, we consider a simplified para- 
metric cycle represented by the boundary condition Q 
with 

a[0) = ^ (cos 6* I cos 6*1), 

b{0) = b'{0) = 0, 

c{0) = ^(-cos6i + |cos6l|), 

d(6l) = ^(-sin0 + |sin6l|). (17) 
The graphs of these functions are plotted in Fig. |4] This 
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III VI 
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FIG. 4. Graphs of functions defined by (17 1 for 9 G [0, 2tt] 



system traverses four distinct segments in the parameter 
space of 6, as 9 runs from to 27r. The segments cor- 
respond to four topologies, cf. Fig. [sj Namely, a ring, 
two disjoint lines, a line attached to a ring, and two rings 
connected at a point. 



IV 



00 



C\D O 



FIG. 5. Cycle of topological changes corresponding to func- 
tions (ITtI. 



QUANTUM ANHOLONOMY 



The system wave functions can be written in the form 



Tpi{xi) — «! sinfcxi + Pi cosfcxi, 
i'2[x2) — a.2 sinfca;2 -t- fii coskxi. 



(18) 





FIG. 6. Anholonomy as a result of the long cycle topology 
change for t — 0.1, s = 1. The ratio L1/L2 is chosen as the 
bronze mean (top), silver mean (middle) and golden mean 
(bottom). 



with 



U = 



V = 



f. 1 

sin kLi cos kLi 

V 

/ 1 

- COS kLi sin kLi 

V 



sin kL2 cos kL2 / 



\ 



■ cos kL2 sin kL2 / 



(20) 



The spectra of the system can be calculated from the 
secular equation 



\AU + kBV\ =0. 



(21) 



In this notation, we have 



It turns out that certain choices of s,t and ratio L1/L2 
result in a quantum anholonomy. Below we present the 
results, separately for the long and the short cycle. 



* = 



















a2 









(19) 



Long Cycle 

We first consider the long cycle of parameter change 
represented by ([t]). We performed a numerical calcula- 
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FIG. 7. Anholonomy as a result of the long cycle topology 
change for t = 0.5, s = 0.5. The ratio L1/L2 is chosen as the 
bronze mean (top), silver mean (middle) and golden mean 
(bottom). 



FIG. 8. Anholonomy as a result of the short cycle topology 
change for t = 0.5, s — 0.5. The ratio L1/L2 is chosen as the 
bronze mean (top), silver mean (middle) and golden mean 
(bottom). 



tion of the spectra of the system for six different settings 
of the parameters s,t,Li/L2. We studied two combina- 
tions of s,i, namely 



• t ^ 0.1. 



1, 



• t = 0.5, s = 0.5, 

and for each of them, we considered three different ratios 
L1/L2 (with Li + L2 = 1), namely 

• the golden mean L1/L2 = (1 + y/5)/2, 

• the silver mean L1/L2 = 1 + \/2, 



• the bronze mean L1/L2 = (3 + vT3)/2. 

We explain the numerical result shown in Figure [6] 
The spectra have period 2tt as a function of 6. This is 
because all the Hamiltonians are 2tt periodic by construc- 
tion. On the other hand, however, some wavenumber do 
not have the 27r periodicity. For example, let us keep 
track of the wavenumber (say ko) of the ground state at 
6* = in the uppermost figure. Around 7r/3 < < 27r/3, 



ko increase suddenly and form three crossings with an- 
other eigenstates. Because the graphs consist of two dis- 
joint parts, these crossings are not avoided. The eigen- 
function corresponding to fco is localized within a disjoint 
fragment, and the other three eigenf unctions involving 
the crossings are localized in the other part. Hence, as a 
result of a 27r-cycle of 0, /cq becomes the wavenumber of 
the third excited state. This is an example of anholonomy 
in eigenvalue and eigenspace. Some other wavenumbers, 
e.g., 4-th excited state at 6* = in the leftmost figure, are 
not involved by the anholonomy. By changing the ratio 
L1/L2, we see various patterns of the quantum anholon- 
omy. Also, the anholonomies appear in various choice of 
the parameters, cf. Figure [7] 



Short Cycle 



Let us proceed to the simplified cycle represented by 



(17), in which four distinct topologies are connected in 



a single sequence. The spectra as functions of the angle 
parameter G [0, 27r) have been calculated for the choice 
of parameters s = 0.5, t = 0.5 and L1/L2 attaining the 
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value of the golden mean, the silver mean and the bronze 
mean. The results are depicted in Fig. |8] 

ON THE BOUNDARY CONDITION 




FIG. 9. Finite approximation of the physical meaning of the 
boundary condition Q for 6 belonging to the sectors I and 
II. The strengths of the 5 potentials in the vertices 1,2,3,4 are 
given by the following formulas. Left: vi = V2 = t*") 
1^3 = W4 = Right: ^Q,V2 = ^^,V3 = ^,v^^\. 




FIG. 10. Finite approximation of the physical meaning 
of the boundary condition Q for 6 belonging to the sec- 
tors III and IV. The strengths of the 5 potentials in the 
vertices 1,2,3,4 are given by the following formulas. Left: 




FIG. 11. Finite approximation of the physical meaning of the 
boundary condition Q for 9 belonging to the sectors V and 
VI. The strengths of the 5 potentials in the vertices 1,2,3,4 are 
given by the following formulas. Left: v\ = " , V2 = jzrj, 
V3 = U4 = Right: vi = ^ > > ^ , 

(ta)^ —ta — 'itsad 1—ta ^. 1 — ta — sd 
■V2 = — , V3 = — ^, «4 = ; . 

The boundary condition ([T]) represents a certain sin- 
gular potential in the vertex that has in general a highly 
nontrivial nature. As it has been shown in [31 [5], the 
physical content of a generic boundary condition corre- 
sponds to a set of infinitely strong 6 potentials sited in 



points in an infinitely small web located in the vertex. 
We sketch how the webs look like for the boundary con- 
ditions of our model in Figures [9l [T0| and pT| In all of the 
finite approximations, the parameter e has to be chosen 
small with respect to the wavelength of the particle. For 
e — > 0, the approximation effectively produces the bound- 
ary conditions Q. The result has been obtained using 
the formulas from [3l H] . It illustrates in a simple way the 
relation between the global topology of the system and 
the topology of the internal infinitely small web. 

DISCUSSION 

The creation of the anholonomy in our model can be 
explained in the following way. In regions I, IV, V and 
VI, when the graph consists of a single component, the 
eigenvalues as functions of the system parameter 9 can 
change their values, but they never cross each other, be- 
cause generically there will be avoided crossings. How- 
ever, in regions II and III, the graph is separated into 
two components (two lines or two rings), i.e., the sys- 
tem is in fact made up of two independent entities, both 
having its own energy levels that can independently vary 
with change of the system parameter 9. Therefore, if the 
parametrization of boundary conditions is conveniently 
chosen, the eigenvalues as functions of 9 can cross each 
other, which leads to the anholonomy effect. 
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